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Low-dimensional representations of exact coherent states of the Navier-Stokes 
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We report that many exact invariant solutions of the Navier-Stokes equations for both pipe and 
channel flows are well represented by just few modes of the model of McKeon & Sharma J. FL Mech. 
658, 356 (2010). This model provides modes that act as a basis to decompose the velocity field, 
ordered by their amplitude of response to forcing arising from the interaction between scales. The 
model was originally derived from the Navier-Stokes equations to represent turbulent flows and 
has been used to explain coherent structure and to predict turbulent statistics. This establishes a 
surprising new link between the two distinct approaches to understanding turbulence. 


The problem of finding simple predictive descriptions 
of turbulence has endured since at least the time of 
Reynolds. Recently, two viewpoints have emerged that 
explain structure in turbulence in quite different ways: 
firstly, in terms of invariant solutions of the Navier-Stokes 
equations, which has been used to explain the transition 
to turbulence; secondly, in terms of selective amplifica¬ 
tion or filtering of a superposition of travelling waves. In 
this paper we show that the latter approach efficiently 
captures the structure of these invariant solutions, pro¬ 
viding a new and surprising link between the two distinct 
approaches and supporting the idea that these invariant 
solutions share the same dominant mechanisms as flows 
in the turbulent regime. 

The first viewpoint comes from treating the infinite¬ 
dimensional Navier-Stokes equations that govern turbu¬ 
lence as a nonlinear dynamical system. The programme 
of work arising from this viewpoint has centred on finding 
invariant solutions of the Navier-Stokes equations that 
appear constant in a co-moving frame of reference m, 
and on finding periodic orbits jug. It is hoped that such 
invariant solutions may eventually be used in a weighted 
expansion to compactly describe turbulent flows |7 . 

These invariant solutions arise in pairs at finite am¬ 
plitude via a saddle-node bifurcation at a particular 
Reynolds number. The so-called lower branch (L) so¬ 
lution of each pair denotes a state with lower drag than 
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its corresponding upper branch (U) solution. These solu¬ 
tions are thought to underlie the structure of turbulence 
by concentrating state space trajectories in their vicinity. 
Although the dynamical systems description has been 
most successful at describing transitional flows, it has 
been argued that such solutions are relevant to turbulent 
flow [i m and recent experimental evidence supports the 
view that these solutions continue to be important in tur¬ 
bulent flows and are ultimately responsible for turbulent 
statistics m- 

The second viewpoint is the model of McKeon & 
Sharma which arose from systems and control theory 
HU H2]. This approach treats turbulence as a super¬ 
position of travelling waves, which are attenuated or am¬ 
plified according to their interaction with the mean flow, 
and excited by other travelling waves. In this model, 
the structure and robustness of turbulence comes from 
the interplay between this linear amplification and an 
energy-conserving nonlinear feedback mechanism. The 
resolvent formulation generates an ordered set of basis 
functions by choosing the velocity fields arising from the 
most amplified forcing, the next most, and so forth. The 
model has been used to make predictions about the spa¬ 
tial organisation of turbulent velocity fluctuations m 
and turbulent fluctuation energy spectra mm- The 
resulting modes are travelling waves with phase and am¬ 
plitude that varies spatially. Unlike approaches such as 
Dynamic Mode Decomposition m, Proper Orthogonal 
Decomposition m, or wavelets [18], the model is derived 
from the equations rather than from an existing data 
set. Notably, this viewpoint is entirely in the frequency- 
domain; kinematic descriptions are abandoned in favour 
of a system-level selection of travelling waves. The ori¬ 
gin of these basis functions has a clear physical interpre- 
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tation. The mechanisms are high amplification at the 
critical layer, where the phase velocity equals the flow 
velocity; the lift-up mechanism, where the flow velocity 
fluctuations extract energy using the shear in the mean 
flow; and high amplification for modes with long stream- 
wise wavelength. 

The presence of only one phase velocity in the invari¬ 
ant solutions used here greatly simplifies the problem of 
comparison to the resolvent formulation, in contrast to 
difficulties encountered in the turbulent case [19]. Thus, 
the frequency-domain view of turbulence as a superposi¬ 
tion of interacting travelling waves is well suited to the 
analysis of exact solutions. 

Both the control theory viewpoint and the nonlinear 
dynamics invariant solutions viewpoint bring different 
and important insights, so unifying these distinct ap¬ 
proaches would be an important advance in our under¬ 
standing of turbulence. In this Rapid Communication, 
we show that the exact invariant solutions often are well 
represented by a relatively small number of model modes. 
This shows that the same mechanisms are dominant in 
the invariant solutions as in the model, and therefore as 
in turbulent flows. 

In the following, we project exact invariant solutions 
in pipe and channel flow onto basis functions (modes) 
generated by the model from the mean velocity profile 
of the solutions. We use the notation U B for the bulk 
velocity, R for the pipe radius, h for the channel half¬ 
height, u r for the friction velocity and v for the kinematic 
viscosity. 

The pipe solutions, presented first, were generated 
by continuation using the pseudo-arclength method to 
Re B = 2 U b R/v = 5300 {Re r = u T h/v = 106 - 214) 
from the solutions of [20] using openpipeflow.org. The 
wall-normal resolution was 60 points. These solutions are 
classified into N-class and S-class. The N-class solutions 
have mirror, shift-and-refiect and rotational symmetries, 
with wavy fast streaks and slow streaks arranged to inter¬ 
act with quasi-streamwise vortices. The S-class have only 
shift-and-reflect symmetry, but are otherwise similar in 
structure. Six S-class solutions and ten N-class solutions 
were used, of which four were upper branch and the rest 
lower branch. The N-class upper branch solutions have 
a friction factor close to that of turbulent flow, whereas 
the others are close to laminar flow. 

The channel solutions are from families PI, P3 and P4 
of [21] and were generated using the code channelflow 
[22] . The wall-normal resolution was 81 points. The 
PI (at Re r = 75) and P3 (at Re r = 85) families are 
active in the core of the channel, and approach lami¬ 
nar as Reynolds number increases. There is as yet no 
widely accepted theory for the mechanism that drives 
these solutions. The P4 solutions (at Re r = 85) are 
highly nonlinear with fluctuations localised near the crit¬ 
ical layer. Their sustaining mechanism is well understood 
m E3]. The critical layer for these solutions varies spa¬ 
tially. The PI and P3 lower branch solutions have been 
continued to higher Reynolds number by the pseudo- 


arclength method. For these solutions, the importance 
of the critical layer mechanism becomes clearer at much 
higher Reynolds number [23H25] . 

The systems model from which the basis functions de¬ 
rive is formulated from the Navier-Stokes equations as 
follows. In the following, the three-component velocity 
field is denoted by U (x, t) and the long time-averaged 
velocity field is denoted by Uo(x), with x being a point 
in the flow interior and t being time. The mean veloc¬ 
ity Uo and associated pressure po are assumed known 
a priori. The fluctuations are then u = U — Uo- The 
Navier-Stokes equations can be put in the form 

_ = ^Vp~Uo-Vu-u-VU 0 + ifcr 1 V 2 u + f (1) 

fo = Uo • VUo + Vpo — i?e _1 V 2 U 0 (2) 

f = —u • Vu (3) 

0 = V • u = V • Uo- (4) 

The model formulation proceeds by considering a su¬ 
perposition of fluctuations in an infinite pipe or chan¬ 
nel, of purely harmonic form at temporal frequency cj, 
streamwise wavenumber cq and azimuthal (spanwise) 
wavenumber /3, allowing the first equation (linear in the 
fluctuations) to be considered as harmonic disturbances 
forced by the interaction between other harmonic dis¬ 
turbances. The phase velocity is then c = uj/a. The 
equation for the fluctuations is then of the form 

u(y-a, /?, c) = K„, ftc f (y;a,/3,c) (5) 

where y is the wall-normal distance and the " notation in¬ 
dicates the appropriate complex Fourier coefficient. The 
object of the analysis is the linear operator which 

is known as the resolvent operator. The analysis then 
considers the singular value decomposition of %, 

= '52^m(y,a,l3,c) a m (a,P,c) 4>* m {y,a,fi,c) (6) 

m 

By definition, the left and right singular vectors and 
the singular values obey the orthogonality and order¬ 
ing conditions, (<t> m (y; a, /?, c), 0 m / (y; a, /?, c)) y = 

(VVn (jj’-i fli C) 5 VVn' (2/; ft , c ))y ~ ^m,m' •> &m &m-\-1 • 

The singular values a m are the amplification factors 
from f to u and the left singular vectors are the basis 
functions (response modes) which represent the velocity 
field. The singular values each represent the gain from 
forcing with the associated right singular vector. This 
gain is assumed to rank the importance of a mode pair in 
a flow, and thus induces a natural ordering of the modes. 

Given the particular Uo for each invariant solution this 
therefore results in response modes particular to each 
solution onto which it may be projected, with a m indica¬ 
tive of the importance of each mode. Only modes with 
the appropriate phase velocity need to be considered. To 
the extent that the modes and singular values correctly 
capture the relevant physics of the solution, only a small 
number of modes will be needed. 
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FIG. 1. (Color online) The upper set refer to the pipe 
solutions (all at Res — 5300), the lower set to the channel 
solutions (at a range of Res)- Left: All invariant solutions 
considered in this study, covering a range of solution classes, 
Res , Re T and wavespeeds c. Right, upper set: Fraction of 
energy captured by a projection of m = 1, 5,10 model modes 
per Fourier mode (■,■,■; pipe solutions). Right, lower set: 
m— 1, 2, 5 model mode pairs per Fourier mode ( channel 
solutions). 


FIG. 2. (Color online) The upper set refer to the pipe so¬ 
lutions (all at Res = 5300), the lower set to the channel 
solutions (at a range of Res)- Left: All invariant solutions 
considered in this study. The quantity e is the internal dis¬ 
sipation of fluctuations relative to the laminar flow with the 
same bulk velocity. Right, upper set: Fraction of total inter¬ 
nal dissipation £ due to fluctuations captured by a projection 
of m = 1,5,10 model modes per Fourier mode (■,■, ; pipe 
solutions). Right, lower set: m = 1,2,5 model mode pairs 
per Fourier mode ( channel solutions). 


A set of exact invariant solutions for channel and pipe 
geometries, broadly representative of all known lower and 
upper branch solutions with single c, were projected onto 
the modes given by the model. 

The efficiency of all the projections of the pipe and 
channel solutions are shown in Figure [l] along with de¬ 
tails of the solutions. The response modes for the channel 
come in pairs that have odd and even symmetry about 
the centreline, so the projections are listed using pairs of 
modes, in accordance with this. Plots representative of 
cases of interest for the solution velocity fields projected 
onto the left singular vectors of the model are shown in 
figures [3] to [6] 

From the projections, we find that all the lower-branch 
pipe solutions, and one of the upper branch pipe solu¬ 
tions, are captured very well by one response mode per 
Fourier mode. In this sense, the model predicts the wall 
normal form of the velocity fluctuations. We also find, in 
particular with the P4 channel solution, that the fluctu¬ 
ation energy is typically concentrated around the instan¬ 
taneous critical layer, where the phase velocity equals the 
instantaneous velocity. This mechanism is known to be 
well captured by the model via the average critical layer 
mm- The extent to which the instantaneous critical 
layer deviates from the average critical layer depends on 


the solution in question. 

The other two upper branch pipe solutions require 
more modes to achieve fidelity. The N4U upper branch 
solution is the most poorly represented solution investi¬ 
gated, with only 80% of the fluctuation energy captured 
by the first ten response modes per Fourier mode. We 
do not know why it is relatively so poorly captured, but 
recent projections of the turbulent attractor show that 
this invariant solution is strongly repelling [26] . It is also 
noticeable that its mean velocity profile looks entirely un¬ 
like that of either the turbulent or laminar flow. The P3U 
solution is also relatively poorly captured. Examination 
of this solution shows that it has a relatively fine struc¬ 
ture, with energy at many Fourier modes. Similarly, we 
see that the solutions with more activity at small scales, 
and so higher dissipation, require mode response modes, 
since the leading modes tend to be smoother. 

Calculations of the internal dissipation of the projected 
solutions show that the true dissipation of the solutions 
is quite well captured for most solutions (Figure [2]). The 
best-represented solutions (in the sense of dissipation) 
capture almost all of the dissipation in the first mode. 
These solutions are close to laminar. The least well 
captured solutions are typically very energetic at higher 
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FIG. 3. (Color online) N3L, lower branch solution in a pipe. From left to right: actual solution; projection onto five response 
modes per Fourier mode (containing 98% of the fluctuation energy); projection onto one response mode per Fourier mode 
(containing 95% of the fluctuation energy); mean velocity profile. The red and blue shading indicates streamwise velocity 
fluctuation faster and slower than the mean velocity, respectively (as a fraction of the maximum amplitude streamwise velocity). 
The quiver arrows indicate in-plane velocity. The wall-normal region where the phase velocity is closest to the mean velocity 
is indicated by a dashed green line in the pipe cross-sections and a red dot indicates the phase velocity in the mean velocity 
profile. The lower branch solutions such as this one are close to laminar, as seen from the mean velocity profile. The mean 

velocity profile (—) has superposed the mean flow generated by the projections with 1, 5, 10 singular values (light to dark-) 

and the laminar (—•). 



FIG. 4. (Color online) N2U (a — 1.25), upper branch solution in a pipe. From left to right: actual solution; projection onto five 
response modes per Fourier mode (containing 96% of the fluctuation energy), projection onto one response mode per Fourier 
mode (containing 84% of the fluctuation energy); mean velocity profile. It is interesting to note that due to the flatness of the 
mean velocity profile, the solutions does not possess an average critical layer. 


wavenumbers and so are more dissipative. These solu¬ 
tions are also the least well captured energetically. We 
suspect that these solutions require higher dissipation 
to stabilise them dynamically, meaning more energy is 
scattered into higher and more dissipative spatial modes. 
The lower-order projections are smoother and fail to cap¬ 
ture this finer structure. Further calculations of the skin 
friction for the pipe solution projections (relative to lam¬ 
inar, not shown) show similar results. 

In the following, we pay special attention to the upper 
and lower branches of the P4 solutions; similar results 
are observed in the pipe flow. As shown in [21], the 
mean velocities of P4U and P4L are respectively simi¬ 
lar to von Karman’s turbulent profile and Virk’s profile 
for maximum drag reduction by additive polymers. In 
addition, the normal or ‘active’ turbulent trajectory lies 
within a region close to the P4U solution while the turbu¬ 
lent trajectory occasionally escapes from this region and 
approaches the P4L solution. 

Figure [7] shows the captured velocity components 
(u,v,w), produced nonlinear forcing (u • Vu), and the 
produced uv Reynolds stress using up to 20 most ampli¬ 


fied response mode pairs in both P4 solutions. In both 
solutions, the first few response modes tend to capture 
the streamwise velocity more than the wall-normal and 
spanwise velocities, because the streamwise velocity dom¬ 
inates the velocity fluctuations. In addition, the same 
order of projection captures a relatively smaller portion 
of u in the high-drag solution (P4U) compared to the 
low-drag solution (P4L). This is because the streamwise 
velocity constitutes a larger fraction of the fluctuations 
in P4U than it does in P4L. 

The model reproduces a much smaller portion of the 
nonlinear forcing. For example, the first 20 most ampli¬ 
fied response mode pairs capture less than 60 and 40 per¬ 
cent of u-Vu for the P4L and P4U solutions respectively, 
see figure [7| To some extent, this is expected since most 
of the nonlinear terms are filtered out by the selectively 
high-gain linear mechanisms in the NSE. In this sense, we 
argue that the response modes capture the necessary por¬ 
tion of the nonlinear terms, and therefore, can better our 
understanding of the scaling role of nonlinearity in the 
NSE. The nonlinear forcing produced by the pipe solu¬ 
tion projections was also compared. Similar trends were 
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FIG. 5. (Color online) P4 lower branch solution in a channel. From left to right: actual solution; projection onto five response 
modes pairs per Fourier mode (containing 92% of the fluctuation energy), projection onto one response mode pair per Fourier 
mode (containing 84% of the fluctuation energy); mean velocity profile. 



FIG. 6. (Color online) P4 upper branch solution in a channel. From left to right: actual solution; projection onto five response 
modes pairs per Fourier mode (containing 77% of the fluctuation energy), projection onto one response mode pair per Fourier 
mode (containing 59% of the fluctuation energy); mean velocity profile. 


observed for the dissipation of the projections. The mean 
profiles resulting from various projections are shown in 
figures [3] and [4j 

Figure [7] also shows how the uv Reynolds stress is rep¬ 
resented by the most amplified response mode pairs. The 
captured normal stresses mx, FT, and ww are not shown 
since they follow the same trends as xx, v, and w fluctu¬ 
ations. We see that uv is captured using fewer response 
mode pairs in the P4L solution compared to the P4U 
solution and that fewer response modes are required to 
sustain the mean velocity in the P4L solution. This is 
confirmed in figure [5] (right) where the mean velocity 
induced by the response mode pairs (dashed lines) are 
compared with the mean velocity from the actual solu¬ 
tion. 

We have shown that the velocity fluctuations in fully 
nonlinear exact invariant solutions can be predicted and 
efficiently represented by a model derived to describe 
high -Re wall-bounded turbulence. This supports the idea 
that the same basic mechanisms are present in these in¬ 
variant solutions as in these turbulent flows. Moreover, 
it should be noted that the model formulation is equally 
suited to representing periodic orbits, which it has been 
argued are likely to be more important in the turbulent 

regime mm- 

The methodology studied will greatly help further de¬ 
velopment of the resolvent formulation of wall turbu¬ 
lence, by providing a simplified environment with a single 
phase velocity in which to study the nonlinear interac¬ 
tions within the model. It was shown that, for the more 
complex solutions, even though the nonlinear terms are 


generally not fully captured by the forcing modes, the 
necessary portion of the nonlinear terms from an input- 
output viewpoint is well captured. This implies that a 
relatively small portion of the nonlinear terms can pass 
through the selective filtering action of the high-gain lin¬ 
ear mechanisms in the NSE. This observation can be used 
to distinguish the ‘active’ nonlinear terms from the ‘inac¬ 
tive’ ones and may have significant implications for mod¬ 
eling and control of wall-bounded flows. 

Because of the small number of coefficients involved, 
it is hoped that it will become much cheaper to solve for 
solutions in coefficient space directly, giving low-order ap¬ 
proximate solutions to exact invariant solutions. Thus, 
we hope that low-order approximate coherent structures 
synthesised from the model will be used to provide seeds 
for the search for new invariant solutions that are already 
close to those solutions. It is hoped this will greatly re¬ 
duce the computational cost of such searches and the 
technical difficulty at large flow rates, and is a subject of 
our future work. 
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(a) 



FIG. 7. The percentage captured for velocity components (u,v,w), for nonlinear forcing (u • Vu), and for the uv Reynolds 
stress using increasing numbers of pairs of leading response modes in the P4L (a) and P4U (b) solutions. 
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